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Abstract 

We discuss the sensitivity of neutrino experiments at the luminosity frontier to generic 
hidden sectors containing new (sub)-GeV neutral states. The weak interaction of these states 
with the Standard Model can be efficiently probed through all of the allowed renormalizable 
'portals' (in the Higgs, vector, and neutrino sectors) at fixed target proton beam facilities, 
with complementary sensitivity to colliders. We concentrate on the kinetic-mixing vector 
portal, and show that certain regions of the parameter space for a new U(l)5 gauge sector 
with long-lived sub-GeV mass states decaying to Standard Model leptons are already severely 
constrained by the datasets at LSND, MiniBooNE, and NuMI/MINOS. Furthermore, scenar- 
ios in which portals allow access to stable neutral particles, such as MeV-scale dark matter, 
generally predict that the neutrino beam is accompanied by a 'dark matter beam', observable 
through neutral-current-like interactions in the detector. As a consequence, we show that the 
LSND electron recoil event sample currently provides the most stringent direct constraint 
on MeV-scale dark matter models. 
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1. Introduction 



The scale at which new physics is to be expected is not always straightforward to ascertain 
from empirical data. The generic arguments for uncovering a Higgs (or Higgs-like) sector 
at the weak scale, in tandem with additional structure that alleviates the fine-tuning issues 
of the Standard Model (SM), leads to the pervading notion that all new physics should be 
found at the high-energy frontier. In part, this results from strong existing constraints on 
new physics charged under the SM gauge group, as would generically be the case in the 
context of electroweak symmetry breaking. However, one of the most compelling empirical 
arguments for new physics is the need to explain dark matter (DM), which suggests instead 
the presence of a hidden sector, uncharged (at low energies) under the SM gauge group. 
The mass scale of the states in a hidden sector, even if the dark matter candidate itself is 
a WIMP with a weak-scale mass, is far less constrained. The notion of a hidden sector is 
also a generic feature of supersymmetric models of the weak scale, due to the specifics of 
supersymmetry breaking and mediation. It is also notable that the presence of new physics 
in a hidden (uncharged) sector is quite compatible with the strong hints toward unification 
of the SM gauge groups at a high scale. 

It is worthwhile recalling that, even in the SM, some of the matter fields are uncharged 
under one or more of the color and electroweak gauge groups, and thus the extension to 
include a further sector is not particularly exotic from a theoretical perspective, even beyond 
the empirical motivations noted above. However, it is clear that direct sensitivity to a neutral 
hidden sector will depend in detail on the couplings, which may be associated with higher 
mass states charged under the SM that have been integrated out, as well as possible light 
neutral states present at low energies. Thus it is helpful to ask in which regimes we have the 
best experimental sensitivity. With respect to the mass and lifetime range for hidden sector 
states, the niche with masses in the range of a few MeV to a few GeV, with lifetimes of less 
than 1 second, is quite intriguing as many of the more severe astrophysical and cosmological 
constraints that apply to lighter states are weakened or eliminated, while those from high- 
energy colliders are often inapplicable. Certain classes of hidden sectors may be strongly 
constrained by flavor physics and other precision tests, but it is clear that many generic 
scenarios are viable. Recent hints from various astrophysical anomalies, when interpreted 
in the context of dark matter, also point to this range which suggests that probing hidden 
sectors with light (approximately) flavor- conserving MeV-to-GeV scale states weakly coupled 
to the SM model represents an intriguing possibility. 

Given a hidden sector, which at low energies necessarily only contains states that are 
electrically neutral and presumably neutral with respect to the full SM gauge group, we can 
parameterize the interactions as follows. 



where O denotes SM and new physics (NP) operators of canonical dimensions k and /, and 
A is a cutoff scale presumably at a TeV or above. In general Onp and O^m can transform via 
various Lorentz representations, but as discussed above we will consider the situation where 
^SM gauge invariant under the SM gauge group. The case of marginal n = interactions 




(1) 



k,l,n 



1 



is especially important and will be the focus of attention here. The SM operators of lowest 
dimension are well known and are often called portals PQ[2]- These operators include: 



F^j, Vector portal (dim = 2), 

H^H Higgs portal (dim = 2), (2) 

LH Neutrino portal (dim = 5/2), 



where F^, H and L are the hypercharge field strength, and the Higgs and lepton doublets. 
The operators ([2]) allow a coupling of the SM to (SM neutral) new physics at the renor- 
malizable level, and thus do not presuppose anything about the mass scale for these new 
fields. Note that the existence of neutrino oscillations can be viewed as new physics entering 
through the neutrino portal, e.g. in the form of a coupling to right-handed neutrinos. It is 
then reasonable to ask whether other portals are realized in nature, and whether we have 
the ability to probe them experimentally. 

Opportunities to probe hidden sectors experimentally arise in particular when states 
exist with a mass scale of a GeV or less and a lifetime longer than the mesons which decay 
due to weak interactions. Such light, long-lived hidden sector states can be searched for 
at high-luminosity fixed target experiments where the large collision-number statistics can 
overcome the (assumed) weakness of the interaction, allowing these states to be produced 
in sufficient numbers. A rough comparison between the collider and fixed target reach for 
new light states in neutral hidden sectors can be made as follows. Assuming the interaction 
between the two sectors is mediated by marginal or irrelevant operators of dimension 4 + n 
with > as shown in Eq. ([T]), the production cross section will typically scale as 



where k is a dimensionless coupling constant and A is a large {e.g. TeV) scale which deter- 
mines the intrinsic weak-coupling 'barrier' to probing the portal. While the characteristic 
integrated luminosity for high-energy colliders is of order lO'^^cm"^, the analogue of inte- 
grated luminosity for a fixed target of length Im with 10^^ protons on target (POT) is 
10^^ X lO^^cm"^ X lO^cm ~ 10^''cm~^. With these numbers, one arrives at the following 
comparison between the production rates for neutral GeV-scale states at colliders and fixed 
targets: 

^y£i^iiil££ ^ 10-6 X (^] " ~ 10-i2+6n^ 

Atarget J 

where in the last relation we have assumed that the attainable collider energies are three 
orders of magnitude larger than s^/"^ for fixed target experiments, E^/ Et = Ed ^J^2/nhpE\^ ~ 
10^ for Ec = 14 TeV and E\sh = 100 GeV. For low n = 0, 1 the production rate at fixed 
targets is clearly advantageous, and even ai n = 2 can be comparable to colliders. However, 
this does not directly imply that the fixed target signal would necessarily be larger, due to 
the geometric acceptance criteria of fixed target detectors. In contrast, if the collider signal 
consists primarily of missing energy, detection may in turn be difficult due to the size of the 
SM background. In essence, (jl]) demonstrates that for n < 2 fixed targets may very well be 
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ahead of colliders in sensitivity to new GeV-scale hidden sector physics, although the precise 
comparison clearly depends on the specific model under consideration. 

Our primary observation here is that the existing and ongoing development of experimen- 
tal infrastructure for probing neutrino physics through beams generated by high-intensity 
proton sources directed on fixed targets is also an ideal setting within which to probe all 
the SM portals, and couplings to a hidden sector discussed above. As we will describe in 
more detail below, past and existing experiments ranging in beam energy from LSND to 
MiniBooNE and MINOS, and those in development such as T2K and NOvA, have strong 
sensitivity to portal couplings. Of course, the idea to search for exotic (un)stable states 
using fixed targets is far from new and in the past several experimental studies have been 
performed targeting different types of exotics [3] . In particular, we note the existing experi- 
mental searches for axions [1], heavy neutrino states [5], long-lived gluinos [6], unstable light 
particles in general [7], and light particles with neutral-current-type interactions [8]. Theo- 
retical motivation for these searches often originated from several unrelated considerations, 
see e.g. [HI UHl [HI [El Il3] for an incomplete list. We stress, however, that the majority 
of these searches were performed in the 80s and 90s with typical POT~ 10^^, which allows 
for significant improvements with modern facilities where the proton exposure can be above 
POT~ 10^^. Moreover, leaving particular theoretical scenarios aside, the description of in- 
teractions in terms of SM portals allows the GeV-scale phenomenology of a hidden sector to 
be systematized. 

In recent years, the scope for particle physics interpretations of dark matter, observed 
through its gravitational impact on cosmology and astrophysics over many time and distance 
scales, has widened considerably. Among numerous dark matter scenarios, several interesting 
possibilities have been proposed that deviate from the most commonly assumed candidate 
of an electroweak-scale WIMP coupled to the SM via higher-dimensional contact operators. 
Notably, models of MeV-scale dark matter coupled to the SM via a sub-GeV mediator 
[m [ISl [ISl [Hi have been proposed in an attempt to explain the 511 keV emission observed 
from the galactic center with annihilating dark matter fl8J. Another suggestion invoked 
WIMPs at the GeV scale |I9] to explain the results of the DAMA experiment [20]. Most 
recently, it was suggested that (sub-) GeV scale mediators between the dark and visible 
sectors could enhance the annihilation cross section of electroweak scale WIMPs to leptons 
[211 [22], in tantalizing agreement with the positron excess observed by PAMELA [23j. While 
many astrophysical anomalies may find alternative explanations, the variety of new effects 
that may occur in dark matter phenomenology due to the presence of light mediators or dark 
matter particles has led to a renewed interest in their low and intermediate energy particle 
physics manifestations [151 [211 [221 [23 [23 [2H]- However, the sensitivity of these tests often 
decreases when the signal contains long-lived particles that escape the detector. In such 
cases, fixed target experiments with detectors lOm-lkm from the target, as in modern long- 
baseline neutrino experiments, can provide complementary sensitivity. A rather striking 
consequence of models with light dark matter is the production of a high intensity 'dark 
matter beam', generated as dark matter particles are pair-produced as a result of the proton- 
target interactions and boosted along the proton beam direction. The scattering of WIMPs 
of this type in the (near-) detector would generate an additional contribution to neutral- 
current type events (see, e.g. [S]). This prediction implies that a direct search for MeV-scale 
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stable WIMP particles is possible at experiments at the luminosity frontier. 

This paper aims to demonstrate that the physics reach of fixed target experiments, e.g. 
in the conventional neutrino beam plus near-detector set up, surpasses the reach of previous 
experimental probes and covers a wider range of possibilities for new physics coupled via the 
vector and Higgs portals. We consider the kinetic-mixing vector portal in detail and estimate 
the number of decays of new GeV and sub-GeV-scale states in past and existing detectors 
such as LSND, MiniBooNE and NuMI/MINOS. In particular, we show that certain parts 
of the parameter space, inaccessible to collider experiments, are already strongly disfavored 
or ruled out by fixed target experiments as a by-product of the neutrino physics analysis. 
Moreover, we demonstrate that DM beam scattering inside the LSND detector places strong 
restrictions on MeV-scale dark matter models, superior to all pre-existing limits. 

The remainder of this paper is organized as follows. In the next section we provide 
some generic details on the coupling of GeV-scale neutral hidden sector states via the most 
important portals. We then proceed to analyze in detail the signatures of long-lived vectors 
and Higgs' bosons from the secluded U(l)5 sector in Section 3. Section 4 contains the 
calculation of the neutral-current-like signal of an MeV-scale DM beam scattering inside 
the detector, while Section 5 provides a discussion of possible future improvements in fixed 
target probes of hidden sectors. 

2. Standard Model portals and GeV-scale states 

Focussing on the portals in ([2]) allows a classification of the interactions of new GeV-scale 
states with the SM. In this section, we discuss the various cases and note regions of param- 
eter space leading to relatively long lifetimes, as appropriate to exploration at fixed target 
facilities. 

2.1 Vector portal 

The simplest way in which the SM may interact with a new abelian gauge sector was pointed 
out long ago ^29j . 

Ant = I v;,F'^^ (5) 

where F^^ and V^v are the SM electromagnetic and secluded U(l)s' field strengths. At higher 
energies, F^^ should be replaced with the hypercharge field strength, but we will ignore the 
coupling to the Z in the energy range considered here. After the spontaneous breaking of 
the secluded 11(1)5, the mixing parameter k appears in front of the coupling of to the 
electromagnetic current, 

/:int = Z.^. + "^hX + ■■■■ (6) 

We have also retained the interaction of the secluded Higgs' {h') and V-bosons, which will 
be important for signatures involving h'. Although the appeal of this model is in its ex- 
treme simplicity, there are nevertheless a considerable range of intricate phenomenological 
consequences (see, for example, [25], [26]). Such models have recently been under scrutiny 
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in situations where the WIMP dark matter candidate is also charged under U(l)5 and the 
hght vector leads to a Sommerfeld-enhanced annihilation rate. 

There are two distinct regimes that allow for a long-lived particle in the U(l)5 sector. 
The first corresponds to k < 10"'^ — 10~^, when the sub-GeV vectors may have a decay 
length CT in excess of several meters. In this case the secluded Higgs' must be heavier than 
the vector, rrih' > my, allowing for a prompt Higgs' decay. Another distinct possibility is 
when K, ~ 10"^ — 10^^ and rriv > rrih'. This is the case of a fast decaying vector, but a 
long-lived Higgs' with a decay width of second order in the mixing parameter k [25l [16] . 
Both possibilities are ideal candidates for a fixed target search, and will be analyzed in some 
detail in the next section. 



2.2 Higgs portal 

Interactions through the Higgs portal, namely couplings to [H [2] are distinct in the 

sense that the Higgs sector of the SM has not been observed, and thus its detailed features 
are unknown. Thus, the Higgs portal is simply one parameterization of an extended Higgs 
sector, and we can take the new state to be a scalar singlet S, that may couple to the Higgs 
portal via dimension 3 and 4 operators, 

Ant = iH^H){XS^ + AS) = hviXS^ + AS) + --- , (7) 

where after the spontaneous breaking of the electroweak symmetry, = (0, {v + h)/\/2), 
and only the interaction terms with the Higgs boson are retained. When A = the model has 
an additional Z2 symmetry so that the particles S are stable and are a viable DM candidate 
The light mass scale for 5* has implications for the natural size of A. Typically, one 



would require Af ^ < m| in order to avoid finely-tuned cancellations. This is not necessarily 
the case for the two Higgs-doublet model with a singlet at large tan/3, as the coupling of 
S"^ to H^Hd does not have to be very small from naturalness arguments alone [32]. Keeping 
to the minimal model, and treating A and A as free parameters, we can integrate out the 
physical Higgs field as we are interested in low energy phenomenology. The result can be 
written in the following form, 

_ (h) Ag^ + AS 

where Cg^ = ^j^fff + ... is the familiar (g^-dependent) SM operator describing the 
interaction of an off-shell Higgs of momentum q with the SM fermions, gluons and photons. 
It is important to keep in mind that unlike the vector case where flavor physics does not 
impose very restrictive constraints on the couplings, the Higgs portal can be constrained via 
rare B and K decays [331 ESI [IS] . 

Long lifetimes for S particles can arise either by symmetry or by kinematic accident. In 
the case of an approximate Z2 symmetry, A <^ v, the smallness of A ensures the longevity of 
S. The lifetime of S is given by the lifetime of the Standard Model Higgs boson at rrih = 
rescaled by the small mixing angle n' = Av/m\: 

Ts = {K')'Tf,{mh = ms). (9) 
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In order to achieve a long lifetime (e.g. in excess of 10~^ s), one has to require a very small 
value for A or to ensure that is below the di-muon threshold. In this case the decay 
to electrons or 7's is suppressed either by the small electron Yukawa coupling or by a loop 
factor. For ms below 100 MeV, the decay lifetime ts is given by, 

100 MeV 

CTs ^ 1 cm X — — . (10) 

Notice that the production of light S states could still be efficient because of the rather 
sizable coupling to nucleons, 

300 MeV/t' ^ - 
CsN ^ X SNN, 11 

V 

so that the coupling constants participating in the production (300 MeV/f ) and decay [me/v) 
can differ by almost three orders of magnitude. 

The new sector coupled to the Higgs portal need not be purely bosonic. Indeed, one could 
use secluded fermions x attached to S scalars [32] , in which case the low-energy Lagrangian 
would also include XX^sm operators. 



2.3 Neutrino portal 

The neutrino portal allows for the coupling of a set of singlet fermions A^^ to the LH compos- 
ite fermionic operator of the SM, and results of course in the neutrino Yukawa interaction, 

C^t = VrjUHNj + ■ ■ ■ , (12) 

where i, j are flavor indices. Together with the allowed Majorana mass terms for the right- 
handed states Nj, these interactions can lead to the generation of observable neutrino masses 
and neutrino oscillation phenomenology, and therefore are well motivated. The flavor struc- 
ture, which is an important aspect of the neutrino portal, is of course the primary means via 
which this portal is probed. Details of neutrino oscillations allow information to be deduced 
about the existence of this coupling even though actual detection of the right-handed states 
is clearly out of reach if these states are heavy, as is motivated in the context of the seesaw 
mechanism for neutrino mass. 

Apart from generating the light neutrino masses, the heavy states Nj can participate more 
directly in SM phenomenology, via the induced mixing angle Or. Given that the right-handed 
states are weakly interacting, they can be produced efficiently only in the weak decays of 
mesons. The typical (naturalness-inspired) mixing angle suggested by the see-saw relation 
is 9'^ ~ m^/Mfi. For Mr ~ 0{GeV) this is still a rather small mixing angle, and detecting 
it would represent a major challenge. Nonetheless, even searches for Nj states coupled via 
larger values of 6 are well-motivated phenomenologically, and have been discussed in the 
past on many occasions [T3] . 



2.4 Higher-dimensional portals 

Models where higher-dimensional operators mediate the connection to a GeV-scale neutral 
sector are ubiquitous. At dimension 5, significant attention has been devoted to axion-like 
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particles, 



■int 



fa 



■^7^75^ H , 



(13) 



where ip represents a generic SM fermion field and the mass of the pseudoscalar a is protected 
against radiative corrections from (fT3|l by the shift symmetry. A departure from the mass- 
coupling relation dictated by the QCD axial anomaly, ~ m^/^//f , towards larger values 
of rria allows the avoidance of various prohibitive astrophysical and cosmological constraints. 
This is a model where fixed target experiments could hold a significant advantage over 
colliders, especially in the kinematic range rria < ^m^ where production is dictated by sizable 
nucleon axial- vector currents, and the decay lifetimes are prolonged by chirality suppression, 
r„ ~ mlma/fl 

Moving up in dimension to dimension six, we encounter a multitude of operators con- 
necting vector and axial-vector currents of the SM with singlet fermions or bosons from a 
hidden sector. A prominent example of such scenarios is the coupling of a light neutralino 
LSP to SM currents. The signal of light neutralino decay, induced either by R-parity viola- 
tion [12] or by a cascade to even lighter super-partners (such as the gravitino) can also be 
studied with fixed target experiments. We will not pursue this subject further in this paper, 
noting that the calculations will in general parallel the pair-production of S-particles from 
the effective Lagrangian ([H]). 

It is also important to bear in mind that there is no intrinsic constraint that the hidden 
sector be weakly interacting. General arguments concerning a UV completion would then 
suggest the presence of nonabelian gauge fields, and thus the natural portal could again arise 
at canonical dimension 6. However, the presence of strong interactions in both visible (i.e. 
QCD) and hidden sectors would allow for the possibility of dangerously irrelevant operators 
that could produce lower dimensional portals in the infrared, i.e. a 'pion portal'. We will 
not discuss such possibilities in detail, but will comment briefly on an example of this type 
momentarily. One of the commonly discussed scenarios of a strongly interacting hidden 
sector is the 'hidden valley' [31], where for example couplings could arise via a pure-gauge 
portal [35] . 




with G the gluon field strength and G' the corresponding field strength in the hidden 
sector. If this is the leading dimension contribution, i.e. the portal starts at dimension 
8, then according to (jl]) colliders would have a clear advantage over fixed target experi- 
ments in probing these sectors. However, in a somewhat more complicated QCD-like se- 
cluded sector, there is the possibility for a dimension 6 operator, connecting quark currents, 
^6 ^(^L(H)7M?i(-R))(^L(_R)7M^L(H))- Such an operator would then translate into the interactions 
of 'hidden pions' with the axial-vector SM current, and so would be essentially identical to 
( |T3|) . In this case, a long-lived 'hidden pion' would be a good candidate for fixed target 
searches. 





(14) 
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3. Fixed target signatures of hidden sectors 



We begin this section with a discussion of production and detection of hidden sector particles 
at proton fixed target experiments. In general there are two mechanisms by which hidden 
sector states Y will be produced: 1) direct production p + A ^ Y + X, and 2) production 
of secondary hadrons p + A ^ H + X followed by the decay H ^ Y + X. In general the 
number of candidate events can be written as 

A^events ^ Ny X Pdet, (15) 

where Ny is the total number of Y particles produced in the experiment and Pjet is the 
probability that a Y particle intersects and decays within the detector. We discuss each of 
these in turn. 

(i) Direct partonic production: For direct production, the number of particles produced 
is simply the cross section times the integrated luminosity: 

Ny = ax iVpoT (^^] , (16) 



where A^pox is the number of protons on target, Na is Avogadro's number, p and A are the 
density and mass number of the material comprising the target, and L is the length of the 
target. For direct production there are essentially two regimes of interest. When the kinetic 
energy of the protons in the beam is large, such that ^/s ^ few GeV, then the production 
cross sections can be computed using perturbative QCD. On the other hand, for low energy 
beams a more appropriate description of scattering is based on the baryonic and mesonic 
degrees of freedom. For order-of-magnitude estimates of the number of events, the subtleties 
of hadronic model-dependence are inessential. However, precise predictions generally require 
a parameterization of the production cross section for a given process, based ideally on data 
from dedicated experiments, as is done for neutrino oscillation studies, as well as on detailed 
simulation of the beam-target interaction and subsequent particle propagation. 

(ii) Secondary hadronic decays: For Y particles arising from the decays of secondary 
hadrons H, we can write 

Ny^NffX BlH^VX- (17) 

To obtain the number of hadrons Nh produced at an experiment, one can calculate the 
production cross section times luminosity directly. However we will often find it convenient 
and perhaps more reliable to base our estimate of Nh on known experimental quantities, 
such as the measured neutrino flux. 

The probabihty Pdet that a produced Y particle reaches and decays or scatters in the 
detector involves a convolution of the angular and momentum distribution of the produced 
Y particle with the decay probability, 

-Pdecay = exp ( —) - exp ( —) , (18) 
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where di{d2) is the distance from the production point to the entry (exit) point of the 
detector, r is the hfetime of Y, and v (7) is the velocity (boost) of the particle. Clearly, 
the convolution is necessary because the distances di 2 depend on the outgoing direction of 
the Y particle (as well as the detector/beamline geometry), and the velocity and boost are 
determined by the Y particle momentum. The integration is restricted to angles subtending 
the detector. The angular and momentum distributions can be calculated in principle, but 
for low energy proton beams this is not always practical due to the complications of the 
strong interaction. Again in this case, it is useful to use a combination of data and theory to 
obtain a parameterization of the distributions, and base predictions on this parameterization. 

To obtain a simple order-of-magnitude estimate of the number of candidate events, which 
we will be content to do for most of the examples in this paper, we can make the simplifying 
assumption that the production of Y particles is isotropic in the c.o.m. frame with a char- 
acteristic c.o.m. momentum Pc.o.m. depending on the details of the process. This leads to a 
simple expression for the probability Pdet- 



where dQc.o.m. ~ 7o'^^iab is the solid angle subtended by the detector in the c.o.m. frame, 
with 7o the boost used to transform between frames, and Az is a typical distance traversed 
in the detector. In the second line above we have made the approximation that the decay 
length is long compared to the distance from the target to the detector. The approximate 
formula f|T9|) also requires a characteristic velocity and boost of the produced Y particles, 
which can be inferred from the kinematics of the decay or production process. 

Taken together, Eqs. fllStllQp are simple formulae useful for obtaining parametric es- 
timates of event rates. In several examples in this paper we have also performed a more 
accurate numerical calculation of the event rates using a Monte Carlo simulation. More de- 
tails on these simulations are given below. We now move on to signatures of specific models, 
beginning with the U(l)s gauge and Higgs' bosons. 

3.1 Sensitivity to a secluded U(l)5 sector 

We will focus our attention on the vector portal and perform a detailed study of the sensitivity 
to such a sector arising from fixed target proton beams. Existing and planned neutrino beam 
facilities provide impressive coverage in the relevant (sub-)GeV energy range: 




(19) 



• mv,h' < 100 MeV - LSND 



• mv,h' < 1 GeV - MiniBooNE 



• mv,h' Z 1 GeV - MINOS, (T2K, NOvA, Project X,. . . ) 
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Depending on the mass scale of the hidden sector portal, the dominant production mech- 
anisms for the minimal U(l)s model are: 



mv,h' ^ 

mv,h' < 400 MeV 



ttO ^ 7V, -fVh' 



(20) 



mv,h' ^ 1 GeV 



p^,uj,(P^Vh',Vn^' 
q + q^V, Vh\ . . . 
q + g^ Vh', qV,... 



iv) 



The final line, relevant for a U(l)5 sector which lies above the hadronic scale, corresponds 
to direct production from partonic interactions. In what follows, we will consider the proton 
beam experiments which provide the best sensitivity to secluded vectors and Higgs' bosons. 

Sensitivity to V: 

If my < ttIt^, the production of V particles in the decays of tt^ is the dominant production 
mechanism. Indeed, the cross section for producing 7r° in proton-nucleus collisions is large, 
and the subsequent decay of 7r° is electromagnetic, so the production of a 1^7 final state, 
although suppressed by /t^ , does not involve the small electromagnetic coupling a. For 
the secluded U(l)s model, kinetic mixing with the photon leads to a decay 7r° — > -^V with 
branching ratio. 



with a very similar expression for r] decay. Another generic advantage of having 7r° production 
is the abundance of photons and electron-positron pairs created as a result. As these photons 
and leptons propagate, they may in turn interact with the target material and if possible 
produce U(l)5' states. Therefore, the intense proton beam has the attractive feature that it 
is also provides a secondary lepton/photon beam. 

Neutral pions vr^ are produced in primary proton-target interactions, in numbers com- 
parable to those of 7r+. In a given fixed target experiment, if the vector resulting from 7r° 
is long-lived, with a decay length on the order of many meters, it may reach the detector 
and decay. Experiments with high luminosity, sensitive to the dilepton decay products, can 
then be used to probe the U(l)5 parameter space. The highest integrated luminosity (10^^ 
POT) was achieved at the Liquid Scintillator Neutrino Detector (LSND) experiment |36j . 
which ran from 1993-1998 at the Los Alamos Neutron Science Center (LANSCE). Thus, as 
we will discuss, LSND puts by far the tightest constraints on the low-mass regions of the 
U(l)5 parameter space. 

According to Eq. (|T71) . we require an estimate of the sample size of neutral pions, N^^o, 
produced in the primary proton-target collisions at LSND. Unlike the charged 7r+ production 
and reaction cross sections, for which dedicated measurements have been made in order to 
gain a precise understanding of neutrino sources, direct experimental information on neu- 
tral 7r° production in proton-nucleus collisions is very scarce, in part because the dominant 
electromagnetic decay products generally do not escape the target. Therefore, instead of 




(21) 
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calculating the 7r° production directly, we will normalize the number of 7r° produced to the 
number of tt"*", which is directly related to the measured neutrino flux. Since the overwhelm- 
ing majority of neutrinos arise from decays- at-rest, we can easily estimate the number of tt"*" 
produced in the lifetime of the LSND experiment: 

(di2iab/47r)i. 

where $1, = 1.3 x 10^^ u cm~^ is the measured neutrino flux at the detector integrated over 
the time of operation, A^et ~ 2.5 x 10^ cm^, and ((ifiiab/4vr)j, 4 x 10~^ is the fractional 
solid angle subtended by the detector relative to the target. This is by far the largest sample 
of charged pions produced in any existing experiment. 

The total number of neutral pions should be within an 0{1) factor of the number of 
7r+ produced at LSND, which can be inferred by comparing the production of 7r° and vr"*" 
in proton-proton and proton- neutron collisions (see e.g. [37]), or by using simple counting 
arguments related to the number of production channels [38] . 

From Eqs. fll5|17|19l) . and using the decay width of V to electron-positron pairs. 




^v^-u = ^^an'myJl -^1 + ^1 (23) 



we can give a parametric estimate of the number of decay events in the regime where the 
travel distance exceeds the distance to the detector: 



^events ~ "^i/^^ot 



an m 



V 



4 / m.^r \ 2 



where Vdet ~ 2 x lO^cm^ is the volume of the LSND detector. One observes that the number 
of events within the LSND detector can potentially be very large. 

We have also performed a more accurate numerical calculation of the event rate. For 
this, we have made the reasonable assumption that the energy and angular distribution, 
dP / {dTj^dVL) , is qualitatively similar to that of the charged pions and have used a parame- 
terization inspired by Ref. |39j. Once the distribution of the outgoing pions is specified, it 
is straightforward to numerically calculate the probability, Pdet, that a U(l)5 vector passes 
through and decays in the detector. To this end, we have written a simple Monte Carlo code 
that produces a neutral pion, decays the pion via tt" •yV, determines if the vector passes 
through the detector, and if so, calculates the probability of decaying within the detector 
according to Eq. f[T5]l . A large sample of events is generated, and the probability Pdet is 
calculated by properly weighting and summing over events in which U(l)s' vectors cross the 
detector. 

The total number of events in which U(l)5 vectors decay within the detector can then 
be written as ^ A 

iVy-decays ~ ^^^^^ X Br^o^^v- X Pdet- (25) 
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Figure 1: Sensitivity of LSND to decays V —>■ e+e^. The light, medium, and dark shaded 
regions indicate more than 10, 1000, and 10^ expected events respectively. The left panel 
shows events due to vectors arising from vr" —>■ decays, while the right panel shows events 
arising from A(1232) NV . 

In Fig. [1] we show a plot of the number of expected dilepton decay events over the n — my 
parameter space, which can be quite significant. The plot shows the result of the full Monte 
Carlo simulation, but the approximate estimate in Eq. (!24l) actually produces results which 
are in agreement with the simulation at the percent level for long-lived vectors. Given the 
large size of the expected event sample, it appears that much of the sensitivity range in Fig. [1] 
can be translated to an exclusion region since electron events from V — > e^e~ decays should 
automatically be included in the recorded LSND dataset, and in numbers inconsistent with 
the well-measured backgrounds from neutrinos and other sources. However, it is also worth 
noting that there are additional kinematic distinctions. In particular, the higher energy 
fraction of resulting pairs would have similarities with the decay-in-flight charged-current 
electron neutrino events [36], for which the background is small, of 0(10). Thus we see that 
the LSND experiment rules out (or is at least sensitive to) a wide range of masses and mixing 
angles, my < m.^ and k ~ 10^^ — 10~^. 

In addition to production through decays of neutral pions, vectors can also be produced 
in the decays of the A(1232) resonance, A NV, which has a branching of 

BrA^w ^ BrA^iV7 xkUI- , (26) 

where 'Bt/^^n^ ~ 0.005. The total number of A particles should be comparable to the total 
number of pions, as we may anticipate that a significant fraction of pions originate from the 
decays of A. Again we use Eqs. fll5lll7|19p to estimate the number of candidate events with 
a vector decaying in the detector at LSND, which is shown in the right panel of Fig. [H This 
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production process extends the sensitivity range to my ~ 300 MeV. 

The MiniBooNE experiment |10] at Fermilab has also accumulated a significant number 
of protons on target, of order 10^^, and involves a more energetic 8.9 GeV proton beam, 
but a larger 540m path to the detector. We have investigated the event rates of vectors 
originating from 7r° and t] meson decays at MiniBooNE and have found the sensitivity to 
be considerably weaker than for LSND due to the lower luminosity and the longer travel 
distance to the detector, which requires smaller k, for the vectors to have a sufficiently long 
decay distance. These two features also affiict the other neutrino experiments in comparison 
to LSND as concerns their sensitivity to V decays. For example, the MINOS experiment 
|42] also has ~ 10^^ POT, with the near detector positioned at a distance of 965m from the 
target. In contrast, the Higgs' can have a parametrically longer lifetime, and thus higher 
energy neutrino beams allow sensitivity to a higher mass range through Higgs' decays, a 
subject to which we turn now. 

Sensitivity to Higgs': 

We now address the production of long-lived Higgs' particles at larger values of k, where 
the dominant decay h' II occurs at order F/j ~ x (loop factor)^. This requires the 
kinematic regime where nih' < my- Further details of the decay channels which determine 
the Higgs' lifetime in different parts of the parameter space are given in [26] , and illustrated 
pictorially in the left panel of Fig. [2l Here we note that specifically for my 3> mh' ^ 2m f, 
the lifetime of the h' is 

[Note that this corrects a typo in Eq. (14) of [26] regarding the dependence on rrih'.] Given 
the boost in the decay, this allows for a decay length of (9(100m) even for n ~ 10~^. At 
first sight a value of k as large as 10~^ can be problematic due to various independent 
constraints, such as the size of 1-loop contributions to (7 — 2 of the muon and the BaBar 
search for exotic resonances in T(3S') decays [4lj. However, it is premature at this point 
to declare that k is restricted to lie at or below 10~^. The measurement of the muon 
anomalous magnetic moment famously disagrees with the SM theoretical prediction at 3a, 
and the existence of a vector with a GeV-scale mass can actually "help" to resolve the 
discrepancy [25]. Furthermore, the BaBar analysis has thus far targeted the decays of T(3S') 
to metastable pseudoscalar particles, and the analysis makes certain explicit assumptions 
that are inconsistent with the U(l)5 modelQ Therefore, a dedicated search for the 11(1)5 
sector with the full (T(4S')) dataset is well motivated [26], [271 EH] . Finally, we note that any 
potential signal involving h' bosons will depend on the secluded gauge coupling a', so that 
an increase in a' may compensate smaller values of k' and lead to sensitivities similar to 
those we present below. Here we fix the U(l)s coupling to be a' = a. 

^In fact, as the leptonic widths of these resonances are three orders of magnitude smaller than the beam 
energy spread, the production of particles from the U(l)s sector is not enhanced by any of the T resonances 
at B-factories, and it is the non-resonant production of Vj and Vh' states that dominates, as e.g. the 
T — > 7I/ decay is strictly forbidden [26l [28] . 
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With these issues in mind, we proceed to consider the sensitivity to Higgs' production 
and decays for k of order 10~^ — 10~^. As outhned above, for masses in the GeV range the 
production of Higgs' particles from the U(l)5 sector may proceed in a number of different 
ways: via the 3-body decays of pseudoscalars, 7r° •yVh'; via the two-body decays of vector 
resonances, p° Vh', etc; and via direct production at the parton level, g + g — > Vh'. For 
completeness, we include the details of the branching fractions and production cross sections 
for these three classes. For three-body decays of 7r°, we have 



Bi^o^^vh' = - dq — — X Br^o^^y.(g ), (2^ 

^ J{mv+mh)^ [Q ~ ^v) 

where Br^o^^y. (g^) is given by Eq. fl2T]) with the replacement niy — > g^, and 



Ty.^Vkiq') = -a'Vq'f ( ^ ) , (29) 



12 q 

where we have defined f{x,y) = a/ A(1, x, y)(A(l, x, y) + 12i;) in terms of the famihar kine- 
matic function A(p, g, r) = + g^ + — 2pg — 2pr — 2gr. For two-body decays of vector 
resonances such as p°, we have 



-2 



We also note that the production of vector mesons in the primary proton-target collisions is 
generally smaller than 7r+ production by a factor of order 100 (see e.g. |37j), which we have 
accounted for in our estimate of the event rate. 

For more energetic proton beams, such as the 120 GeV NuMI source at MINOS [42j . 
direct Higgs' production at the parton level can also be significant. The parton level cross- 
section for gg Vh' takes the form, 

= (l - ^) " / ^) . (31) 

where s ~ sxiX2 is the center-of-mass energy of the two partons with momentum fractions 
Xi and X2- The total cross-section for a proton colliding with a nucleus of atomic number A 
is a{p + A) = Za{p + p) + {A — Z)a{p + ra), where individual hadronic cross sections require 
convolution with the parton distribution functions (pdfs), 

a{p{Pi) + N{P2) ^V + h' + X) (32) 

= / dxj dx2y2ff{xi)ff\x2)-a{qf{xiPi) + qf{x2P2)-^V + h'), 

Jo Jo j: 

in which denotes either a proton or a neutron. We have used the CTEQ6.6M pdfs |13] 
in the evaluation of the cross sections. The factorization scale was chosen roughly at the 
center-of-mass energy for the NuMI beam/target system. We find that the cross section for 
direct production varies from O.Olpb - 2pb as the Higgs' mass drops from 2 GeV down to a 
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Figure 2: On the left we show the Higgs' decay length cr for k = 10~^ and a' = a. The 
regions are cr < Im (white), 1 - lO^m (light), 10^ — lO^m (medium), and > lO^m (dark). On 
the right, we show the sensitivity of LSND (vr^ '-/Vh' - dark), MiniBooNE {p,uj h'V - 
light), MiniBooNE {p^oo h'V, muon events only - medium), and MINOS (QCD Vh', 
muon events only - medium), indicating more than 10 events expected in the detector. 

few hundred MeV. We should also note that direct production of vectors is also possible via 
the resonant process qq V, but given the sensitivity reach in k the decay length of the 
vector is always too short to be relevant for the experimental setup at MINOS, where the 
near detector - approximately a 7m long cylinder with radius Im - is located 965m from the 
graphite target. 

Using these production processes and the details of the Higgs' lifetime, we have explored 
the sensitivity of a number of neutrino experiments to Higgs' decays in the detector. The 
results are shown in the right panel of Fig. [2] for k = 10~^. These estimates have been 
obtained using the approximate formulae in Eqs. (I15m9p . We have also checked these results 
for the case of LSND and NuMI/MINOS by running a Monte Carlo code that incorporates 
the required branchings/production cross section and kinematics, the decay length of the 
Higgs', and the geometric acceptance of the detectors, obtaining very similar sensitivities to 
those depicted in Fig. [21 Given the distance to the (near) detectors, the proton sources at 
LSND, MiniBooNE and NuMi/MINOS do not have enough luminosity to probe much lower 
values of k, via the Higgs' channel. However, we observe an important complementarity with 
searches for vectors, in that the more energetic proton sources allow this channel to probe 
masses up to the few-GeV range from dimuon events in the MINOS near-detector. 

One final process we will briefly discuss is h' — V conversion in the detector. If the decay 
length is longer than the distance to the detector, it is possible for the h' to scatter off a 
nucleus or electron and convert to a vector, which instantly decays to a pair of leptons. We 
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have estimated this signal and there are indeed a large number of events expected at LSND 
for K = 10~^. However, the sensitive kinematic range is limited to lighter vectors and Higgs' 
bosons, essentially as the Higgs' must have sufficient kinetic energy to upscatter into the 
heavier vector (recall my > rrih' in order for the Higgs' to be long-lived). We conclude that 
LSND is sensitive to Higgs' conversion, but in a region of parameter space already covered 
by h' decays. Moving to smaller values of k, the rate of conversion in the detector is too 
small to expect a significant number of events. 

Thus far we have considered the minimal vector portal scenario, which is justified because 
V and h' comprise a primary set of states forming a module in any more complex hidden 
sector, and thus lead to the most generic predictions for experimental sensitivity. However, 
it is also worth considering which extensions might lead to a sensitivity reach in parts of the 
parameter space, e.g. k, ~ 10~^ — 10^^, which fall somewhere in between the natural range 
for V and h' decays. One obvious possibility would be a hidden sector with at least two 
new gauge bosons Vi and V2, e.g. as part of some nonabelian gauge group, or simply two 
distinct U(l)5' sectors. If the heavier vector state were produced on-shell and subsequently 
decayed V2 Vih', with the Higgs' long-lived, this would provide an additional window into 
the parts of the parameter that are inaccessible in the minimal model. 

Probing short-lived U(l)s particles using a beam: 

Another means of producing long-lived vector states arises via a kaon beam, that may be 
obtained in many fixed target experiments. Previously, the production of V bosons in the 
decays of K"^ was calculated in [25]. Here we note that the rare decays of Kl can also be 
used to constrain the U(l)5 sector; these decays occur well away from the target and may 
therefore have very low backgrounds from other sources. This opens up the possibility of 
probing k, ~ (9(10~^ — 10~^) directly in the decays of Kl. For example, the branching ratio 
for Kl V'-f — > l^l~'j is given by 



This rare decay can be probed by searching for resonance structures in Kl /i^yU^7, which 
has a SM branching ratio of 3.6 x 10"'', measured recently by the KTeV collaboration ^3] 
with almost 10^ events. The search for additional resonant structures among these events 
can also test the vector candidate [25] interpretation for the "HyperCP resonance" [45]. The 
KTeV collaboration has also made a precise measurement of the 7r° e^e~ decay with the 
branching (6.44 ± 0.25 ± 0.22) x lO"*^ [46], and about 800 events coming from the 7r° decay 
produced in Kl 37r°. According to Eq. fl2Tl) . a relatively light vector with ~ 0(10"'') 
should produce roughly the same number of events and therefore can be searched for using 
the same dataset. This would represent a significant gain in sensitivity over B-factories, 
albeit in a rather limited kinematic region. 

Constraints from hidden sector WIMP dark matter. 

One of the primary motivations for anticipating a hidden sector is the need to explain 
dark matter, and while exploration of the portals is formally independent of other features of 




(33) 
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hidden sector physics, it is useful to consider the independent constraints on the parameter 
space that may arise when the dark sector contains WIMP dark matter x- If the dark matter 
candidate is charged under U(l)5', then vector exchange-induced form-factors such as the 
charge radius can lead to significant scattering cross-sections on nuclei, and direct detection 
limits may constrain n < 10~^ even for my ~ 1 GeV. However, even a mild splitting of 
the charged components significantly ameliorates this constraint [17], and so much larger 
kinetic mixing is also perfectly consistent with adding a hidden sector WIMP. A number of 
other astrophysical and cosmological constraints have been considered recently in connection 
with the presence of a dark matter candidate charged under a hidden U(l). We note that 
vectors should decay before BBN, limiting the lifetime to 0{ls). Constraints on the diffuse 
extragalactic X-ray background also impose limits on my/m^ through upscattering of CMB 
photons due to the annihilation xx ~^ W integrated over all halos [IE] . The corresponding 
impact of energy injection during recombination on CMB anisotropics and polarization has 
also been argued to limit my/m^ in these scenarios to be at or above 10~^ [49]. Thus for 
a 100 GeV WIMP charged under the secluded U(l), it seems clear that from a number 
of constraints, the vector mass should be above about 10 MeV. Annihilation of WIMPs 
accumulated in the solar interior into metastable mediators will in addition provide strong 
sensitivity to lifetimes in the millisecond range for metastable U(l)5 particles. 



3.2 Sensitivity to the Higgs and higher-dimensional portals 

In this section we estimate potential sensitivity to light scalars and pseudoscalars coupled to 
the SM via the Higgs portal, and via higher- dimensional operators. Since this is a subject 
which has a rather vast scope, we limit the discussion here to the production of relatively 
light 0(100 MeV) scalars and pseudoscalars. 

The production rate for pseudoscalar particles a with masses comparable to the pion 
mass, and with derivative interactions f[T^ to nucleons, can be obtained by rescaling the pion 
production cross section if one is willing to disregard differences in the isospin dependence. 
In this way, we find the total production rate for C(IOO) MeV pseudoscalars to be 

K-Kx&c. 10-1° X X (^^)' , (34) 

and therefore we can normalize the pseudoscalar flux through the detector to the neutrino 
flux from decay- at-rest. The decay rate of pseudoscalars is dominated by the decays to 
electrons (and photons at the loop level) and is given by. 



mliria 



2vr/, 



(35) 



With the optimal choice for cTq (/„ ~ 20 TeV), there can be up to a billion a decays in the 
LSND detector! In the regime where the decay length exceeds the distance to the detector, 
the expected size of the signal for a 100 MeV pseudoscalar is 

AT ^ ^det /lOTeVy 

Na decay ~ X 10^ X , (36) 

PalaCTa \ fa J 
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where we took /3a7a ~ 2. This estimate shows that even a scale as high as 500 TeV can 
be probed via the decay of axion-hke particles. This event rate can also be translated into 
constraints on "pions from hidden valleys" , should they have a similar kinematic mass range. 
More specifically, the effective decay constant fa can be identified with fhn/-^l, where 1/Ag is 
the normalization of the effective operator connecting the axial-vector currents in the visible 
and hidden quark sectors and fhn is the "hidden pion" decay constant. 

Finally, we discuss the sensitivity to scalars S via coupling through the Higgs portal. 
This topic is closely related to hadron physics in the O"*" channel, and as such is shrouded 
in controversy. At LSND energies, S'-particle production cannot occur from the decay of 
hadronic resonances. In order to obtain an estimate, we use a simple-minded factorization 
formula for the process p + N —>■ p + N + S, together with the coupling constant from 
Eq. ([II]): 

as ~ , , X atot =^ Nsr^ 10"^ x {n'f x N^. (37) 



Air 

The decays for ms — 100 MeV are given by the same formula as for pseudoscalars, Eq. fl5^ . 
swapping k'/{2v). The optimal decay range for the scalars, maximizing the LSND 

signal, is (k')^ ~ 10~^, and therefore the number of decay events from scalars coupled to 
the SM via the Higgs portal can maximally reach (9(10''). Although this is an impressive 
level of sensitivity, the rare decays of K mesons [I6] can probe (k')^ ~ few 10^^ which 
suggests rare decays may be a few orders of magnitude more sensitive than fixed target 
probes. MiniBooNE and other more energetic experiments can access an additional channel 
for producing scalars via rj — > tt^S decays. This branching is enhanced because of the small 1] 
width, that decays either electromagnetically or via three pions with a reduced phase space. 
A detailed analysis of this possibility falls outside the scope of this paper. 

Finally, we comment on the pair production of S particles due to the effective Lagrangian 
coupling ([8]), followed by the decay of the scalars to leptons. The production rate in this case 
is generally very small, suppressed by the fourth power of the Higgs mass in the denominator. 
However, one of the interesting modes is SS, which proceeds via a top-W^ loop and 

benefits from the large top Yukawa coupling. For A ~ 0(1), the branching of to a pair of 
scalars can reach 10~®, and lead to up to a thousand leptonic decays of 5* in the MiniBooNE 
detector if the decay length is optimized by a suitable choice of k'. 



4. Probing a dark matter beam 

An intriguing extension of the idea to probe light states in hidden sectors via one of the 
renormalizable portals is to consider the possibility that the dark matter candidate may 
itself be light enough to be produced in fixed target experiments. This would generate a 
"dark matter beam" which may register in the neutrino detector via neutral current-like 
interactions. To demonstrate this point we will calculate the expected signal at LSND in the 
MeV-scale dark matter scenario, mediated by relatively light vector particles. Such models 
were first studied in [HI [I5] and rephrased in the context of kinetically mixed U(l) models 
in [in]. The original motivation for considering dark matter candidates in this mass range 
come from the seemingly inexplicable strong and spatially homogeneous emission of 511 keV 
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photons from the galactic bulge [50]. It was proposed that MeV-scale WIMPs annihilating 
to electron-positron pairs could in principle source this emission [TB] . 

One of the viable models posits scalar dark matter flQ\ charged under the secluded U(l)5 
having kinetic mixing with photons. One considers a regime distinct from that considered 
above with mv > m^, where x stands for the dark matter scalar. The fact that the DM 
particles are scalars ensures that annihilation proceeds in the p-wave, which is rather impor- 
tant as it allows the use of the same mechanism for the annihilation at freeze-out and for the 
galactic annihilation to positrons, which has to be much smaller. The freeze-out requirement 
translates to 



(}^)\(^y^m-^ (38) 

a \ mv J VMeV/ ^ ' 

From now on, for definiteness we take the WIMP mass = 1 MeV in order to be fully 
consistent with the limits on the shape of the 511 keV line [51] and with constraints on the 
accompanying 7-emission [52]. 

An important feature of this construction is the dominance of the invisible width of V 
to a pair of DM particles relative to the visible width to leptons, 

^ = ^ « 1. (39) 

i inv (-^ 

At first sight this should greatly reduce the chances for observing the products of V decays, 
as opposed to models where the decays of V are dominated by the branching to SM particles. 
However, this need not be the case in situations where there is some advantage in having a 
long-lived (or in this case stable) beam of particles that can scatter in the detector. 

In what follows we estimate the signal in the LSND detector from the following sequence 
of events: 

1. p + p-^ X + TC^ 

2. 7r° ^ -fV 

4. x + e^X + e 

The first three processes in this sequence occur inside the target and the beam stop, while 
the last process by which the energetic electron is created, occurs within the detector after 
the dark matter beam has propagated from the target. The model favors a light vector 
mediator, which could therefore be produced in pion decays. Using the same assumptions 
as in the previous section, and neglecting the phase space factors, we arrive at the following 
estimate for the integrated DM flux through the LSND detector: 

~ (40) 

In order to obtain the number of neutral current-like electron scattering events within the 
detector, we calculate the scattering cross section directly: 

daex^ex ci;''*^ 8na^me{l — Ef / E) 
dEf {m\ + 2m,E^Y ' 
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Figure 3: Expected number of neutral current like electron events induced by MeV-scale dark 
matter scatterings in the LSND detector. We show regions with greater than 10 (light), 1000 
(medium), and 10^ (dark) expected events. Regions below the black line correspond to a 
strong coupling regime, with values of a' > Att. 

where Ef and E are the energies of the recoil electron and incident dark matter particles 
respectively. 

If we conservatively set my ~ 10 MeV and a'/t^ ~ 10~^, then with E ~ 100 MeV we find 
the total scattering cross section to be aex^ex ~ 10~'^^cm^. Combining this with the total 
WIMP flux, the number density of electrons in the detector He, and the detector volume, we 
arrive at the following estimate for the total number of neutral current-like events: 



This result is three orders of magnitude above the total number of energetic electrons pro- 
duced by the neutrino beam and observed in the LSND detector. In Fig. [3] we show the 
number of expected events with electron recoil energies above 20 MeV in the K-my parame- 
ter space. We also show the line below which the dark sector is strongly coupled, a' > Air. In 
fact, constraints on the WIMP self-scattering cross section imply a coupling well below the 
perturbative bound, which further restricts the parameter space. Therefore, we can safely 
conclude that this particular model of MeV dark matter is ruled out by the LSND dataset. 
Even more stringent constraints can be obtained if one considers spectral information on the 
recoiling electrons. Since the DM beam energy is not attenuated at all by the target, the 
energy of the recoil electrons can be larger than the main background produced by neutri- 
nos from stopped pions and muons, and the number of background decay-in-fiight events is 
estimated to be of C(10). It is important to investigate whether the estimate (H2|) can be 
extended to all MeV-scale dark matter models, which would essentially rule out all analogous 
particle physics scenarios of this type, at least in the form motivated by the galactic 511 keV 
line. A possible modification of the model that may escape the stringent LSND limits would 



NnC =neX Vdet X (T, 




(42) 
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involve a kinematic constraint so that no V particles could be produced. In this case, more 
energetic experiments such as MiniBooNE should provide valuable constraints on a dark 
matter beam. A detailed analysis of the fixed target constraints on MeV-scale dark matter 
models will be presented elsewhere |53] . 

5. Conclusions 

In this paper we have explored the sensitivity of the already expansive fixed target experimen- 
tal neutrino physics program to a broader range of portals coupling the SM to hidden sectors. 
We have emphasized that the remaining vector and Higgs portals are also well-motivated 
points on which to focus searches for new physics, particularly through their possible connec- 
tions to the physics of dark matter and other hidden sectors such as supersymmetry breaking. 
As elaborated in sections 3 and 4, existing data from LSND, MiniBooNE and NuMI/MINOS 
is already very constraining for a number of scenarios of new physics entering through the 
vector and Higgs portals, and particular scenarios for MeV-scale dark matter, motivated by 
the galactic 511 keV signal, are essentially ruled out by limits imposed by the LSND dataset. 

In this concluding section we will comment briefiy on other fixed target probes and ex- 
pectations for future improvements. We have focussed our attention on the neutrino sources 
with the highest integrated luminosity in the relevant energy range. However, there are 
several other experiments that explore different kinematic regimes. The NuTeV experiment 
at Fermilab had a very energetic 800 GeV proton beam, and actually observed 3 events 
consistent with a long-lived state decaying to two muons, but with asymmetric momenta 
|54] . However, with only 10^^ POT and a 1.4km distance to the target, this would not be 
consistent with the minimal U(l)5 scenarios considered here. For the lower mass range, i.e. 
my ^ 100 MeV, electron beam dumps can also provide a competitive source of constraints, 
as recently discussed in the context of the vector portal in jSH]. In particular, E-137 at 
SLAG [56j can provide a reach for low mass vectors down to small values of n which is 
competitive with the limits discussed here from LSND. In addition, E-141 (SLAG) and E774 
(Fermilab) provide constraints for larger k due to the small travel distance to the detector, 
but only for very low mass vectors [55]. For the heavier mass range, and with longer decay 
distances, it seems that proton sources hold a clear advantage. Therefore there is significant 
complementarity between electron and proton beam probes of hidden sectors. 

While existing facilities already impose significant constraints, the development of long 
baseline neutrino experiments with high luminosity proton sources suggests significant im- 
provements in sensitivity in the near future. The NuMI beam is already in use for the MINOS 
experiment, and will be upgraded for use with NOvA which began breaking ground on the 
site for the far-detector in May 2009. Meanwhile, beam commissioning for T2K began in 
April 2009. A notable feature of T2K in the present context is the possible installation, in 
addition to the near detector and beam monitors located 280m from the target, of a second 
near-detector at a distance of a few km that would allow a combined sensitivity to a larger 
mass range for metastable particles up to a few GeV. In conclusion, further development of 
the experimental programme devised to study the physics of neutrino fiavor can at the same 
time play a significant role in probing the remaining renormalizable portals, and thus will 
have a primary role to play in exploring hidden sectors containing at least some states which 
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are light relative to the weak scale. 
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